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ABSTRACT 

We present simulated results of quasi-periodic flares generated by the inelastic 
collisions of a star bound to a super-massive black hole (SMBH) and its attendant ac- 
cretion disc. We show that the behavior of the quasi-periodicity is affected by the mass 
and spin of the black hole and the orbital elements of the stellar orbit. We also evaluate 
the possibility of extracting useful information on these parameters and verifying the 
character of the Kerr metric from such quasi-periodic signals. Comparisons are made 
with the observed optical outbursts of OJ287, infrared flares from the Galactic center 
and X-ray variabihty in RE J1034-f 396. 

Key words: accretion discs, radiative transfer, galaxies: active, galaxies: individual: 
RE J1034-I-396, BL Lacertae objects: individual: OJ 287, Galaxy: center. 



1 INTRODUCTION 

One of the important developments in our understand- 
ing of galactic nuclei has been the realization that stars 
are intimately involved in the dynamics of accretion discs. 
Direct evidence for this comes from observations of Sgr 
A*, where three-dimensional stellar orbits are determined 
within lO^^cm of the compact rad i o source with no pre- 
ferred orientation llChez et al.1 l2005l : ICenzel fc KarasI l2007l : 
iGillessen et al.ll2009h . These observations lend support to 
the notion that stars can actually be formed from accret- 
ing gas in a disc. This realization leads to consideration of 
the dynamical interaction between stars and discsand the 
observable consequences of this interaction. 

The passage of a star through an accretion disc has 
been the focus of many studies. Stars can influence the 
disc by impo sing a drag and removing its angul ar mo- 
mentum (e.g. IOstr iker"l98^: ' Norman fc Silkl Il983l ). heat- 
ing (e.g. Perrv fc Willia ms 1993,), and depositing mass (e.g. 
lArmitage et al.l Il996l ). On the other hand, the distribu- 
tion of the stellar orbits can be altered by the disc (e.g. 
ISver et al.|[l99ll : lRauchlll995l : iKaras fc SubiH200ll ) . IZentsoval 
showed that each passage through a disc results in 
the appearance of a bri ght hot spot on the disc's surface. 
iPerrv fc WiUiamd (|l993l ) related this passage to a UV bump 
superimposed upo n the continuum sp ectrum emitted by an 
optically thick disc. lZurek et al.] l| 19941 ') showed that star-disc 
collisions can leave "star tails" above the disc and produce 
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broad line emission. iNavakshin et al.] (|2004l ) used the pas- 
sage to explain a large mag nitude and short dura tion X-ray 
flare at the Galactic Centre (|Baganoff et al.ir200ll ). However, 
a bound star will pass many times through a disc, each pas- 
sage leading to a flare, as shown in Fig.[T] The timing of these 
flares can therefore be a measure of the elements of the stel- 
lar orbit and, when the orbit is relativistic, it can provide a 
probe of strong-field gravity. In this paper, we consider what 
information could, in principle, be derived from careful ob- 
servations of a single star in an individual relativistic orbit 
about a black hole and its attendant disc. 

In classical mechanics, when a star orbits a massive 
black hole, the stellar trajectory is a bound ellipse with 
the position of the heavier object at the focus. However, in 
the strong gravitational potential well of black hole, a star's 
trajectory will be more complicated. The three-dimensional 
orientation of each stellar orbit differs due to two general 
relativity effects. First, the peribothon, or closest approach 
to the hole, precesses on the orbital plane. Second, if the 
black hole is spinning, due to the dragging of inertial frames, 
the orbital plane itself will precess about the spin of the 
hole; henceforth, t he orbit becomes even more complex (e.g. 
iMisner et alll 19731 ). 

Not only massive particles but also massless photons' 
trajectories are altered by space-time curvature. Relativis- 
tic ray-tr acing has been tackled both an alytically and nu- 
merically. ICunningham fc BardeenI l| 19731 ) use this method 
to derive the transfer function that allows one to im- 
age objects near black holes. The propagation of pho- 
tons from the accretion disc around black holes has been 
well studied using semi-analytical methods and dire ct nu- 
merical methods by many people (see for example, IPiranI 
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Figure 1. A star bound in the gravitational potential of the black 
hole: it orbits around the black hole, hits the accretion disc, and 
produces flares repeatedly. 
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Combining the stellar trajectory, the star-disc collision 
model, and ray-tracing the photons from the flares, the time 
behavior and flux of the flares caused by star-disc collision 
can be modelled. These signals exhibit quasi-periodicity for 
two reasons: (1) due to the precession of a stellar orbit, in- 
trinsic time separations between two consecutive star-disc 
collisions are different, and (2) because the star hits the ac- 
cretion disc at different (r, 0) coordinates on the equatorial 
plane each time, the light- crossing time for the signals is also 
different. 

This paper is organized as follows. In section 2, we 
give the assumptions and methodology, describe the stellar 
trajectories, construct our collision model, and discuss the 
method of ray-tracing. In section 3, we present simulated 
results for the quasi-periodic signals, showing in detail how 
the quasi-periodicity depends on the parameters of the black 
hole and the star's orbit. In section 4, we discuss the possi- 
bility and method of extracting useful information from the 
signals. In section 5, we relate our model to observations of 
quasi-periodic behaviors in AGN. Further implications are 
briefly discusses in the last section. 



2 MODEL AND METHODOLOGY 

2.1 Assumptions and methodology 

The prime assumption in this paper is that the gas particles 
in the accretion disc and the star orbiting the SMBH can 
be treated as massive point particles. (If there is significant 
pressure in the disc, the influence on what follows will be 
minor.) Consequently, their trajectories are geodesies in the 
space-time of the black hole. The gravity from the star and 
disc are assumed to be inconsequential compared to that of 
the hole, so the vacuum metric may be used to calculate 
these trajectories. 

The accretion disc around the black hole is assumed to 
be thin and Keplerian. The outer radius is set as 20 gravita- 
tional radii (Rg). In this paper we use units with G = c = 1 
for simplicity, which gives Rg = M (the mass of the black 
hole). S = J/M is the angular momentum per unit mass for 
the black hole, and S = S/M is defined as the spin, which 
fall within the range between -1 and 1. 5* = corresponds to 
a spinless Schwarzschild black hole and S = 1 corresponds 



to a maximally rotating Kerr black hole. The disc's inner 
radius is the innermost stable circular orbit (ISCO), which 
is at 9M, 6M, and 1.237M respectively for black holes with 
S = -1,0, and 0.998. 'We take spin S = 0.998 in stead of 
5= 1 as a reasonable upper bound of spin (e.g., iThornd 

Hizi)). 

The collision of the star and the gas transfers part of the 
star's kinetic energy to heat or radiation. This energy per 
orbit lost is small compared with the kinetic energy of the 
star, leaving the motion of the star roughly unaffected. We 

also assume we are well outside the tidal radius: r 3> 
for main sequence stars. 

The trajectory of a star is fully determined by its ini- 
tial position and momentum. Eight initial conditions are re- 
quired: to,ro,9o,4>o,io,ro, ^o, and (po- T he energy E, angular 
momentum L and Carter constant Q (ICarteij[l968l ) can all 
be derived from these initial parameters. The observer is 
defined to be at some inclination 9 to the black-hole spin 
direction looking down on the accretion disc at the flares. A 
fifth-order Runge-Kutta numerical method is used to calcu- 
lated how a particle moves given its position and momen- 
tum. 



2.2 Thin Keplerian disc 

'We use the Boyer-Lindquist metric to describe the SMBH 
space-time: 



+T.d6' -f + 5 



E 

2S'^ Mr sin^e 



A 



sin^ edfj)^ , (1) 



with E = -f 5"^ cos^ 9, A = r^ - 2rM + . 

We assume that the accretion disc is equatorial and 
Keplerian. For such a disc, its gas particles will have an 
angular velocity: 




(2) 



(|Misner et al.ll 19731 ) for a prograde orbit. 

Combing with the condition that the inner product of 
the four velocity of the gas should be -1, we have: 




where 



(3) 



(4) 



These equations describe the four- velocity of the gas parti- 
cles in the disc. 



^ A magnetized neutron star provides an interesting possibility 
as it can present a large effective area close to the horizon of an 
intermediate mass hole. 
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2.3 Inferring stellar orbital elements from 
observations 

We treat stars as massive point particles following geodesies. 
The trajectory of a bound star in general relativity is ap- 
proximately an ellipse, described by the orbital elements: 
(the semi-major axis of the orbit), e (the eccentricity), i (the 
orbital inclination angle with respect to the accretion disc), 
Q (longitude of the ascending node N with respect to the 
projection of the hue of sight on the disc), uj (angular dis- 
tance between N and peribothon), and To (the epoch). The 
three-dimensional orientation of stellar orbits varies due to 
two relativistic effects. Firstly, the peribothon precesses in 
the orbital plane. Secondly, if the black hole is spinning, the 
orbital plane will precess about the spin of the hole due to 
the dragging of the inertial frames. 

Figs. [2] and |3] show the trajectories of stars with ar- 
bitrary initial parameters with black hole spins S=0 and 
S^COOS, respectively. We can see that in the S—O case, the 
star's closest approach to the black hole precesses yet the 
motion remains on the same plane, forming a rosette pat- 
tern. By contrast in the S'^O.OOS case, the star's orbit is no 
longer planar - the orbital plane precesses around the spin of 
the black hole as well. The resulting trajectories look more 
like a ball of tangled yarn, with the amount of extra pre- 
cession depending on the peribothon and the inclination of 
the stellar orbit, and the spin of the hole. For both setups, 
we can see that precession is rapid close to the hole, thus 
making the orbit quite irregular. 

The mass of the hole and the size of the stellar or- 
bit need to satisfy different conditions in order to observe 
general relativistic features of the stellar orbit around the 
SMBH: e.g. produce enough flares in reasonable observing 
time, reveal the eccentricity through the apsidal motion, re- 
veal the spin through the Lense-Thirring effect, have a grav- 
itationally stable stellar orbit, and avoid tidal disruption. 
If these constraints are satisfied, we may be able to derive 
the parameters of the system from the temporal behavior 
of the signals. This provides, in principle, a way to measure 
the spin of the black hole through accurate timing measure- 
ments. 



2.3.1 The size bound 

The ISCO describes the closest stable circular orbit around 
the black hole. For a black hole with mass M and spin 5", 
the semi-major radiu s of the orbit ttr satisfies the equation 
l|Bardeen et al.lll972D : 



35" = 0. 



(5) 



Therefore, flr increases linearly with M for a given S. If we 
choose a maximal value of 5* = 0.998, we get the lower limit 
constraint on the size of a circular equatorial stellar orbit 
for the star not to plunge i nto the event horizo n of the black 
hole. This corresponds to ijMisner et al.l[l973l ) 



Ur > 1.237M 



(6) 



for a circular stellar orbit around a spin 0.998 black hole, 
which in cgs unit is: 



M < 5 X 10' 



/ ar \ 
UOi^cm/ 



Ma 



(7) 



2.3.2 The period bound 

If the period is too large, then useful observations will take 
too long. Therefore, we set the orbital period P to be less 
than 1 year to obtain enough fiares in a reasonable amount 
of time. 

Here the spin term is relatively unimportant when the 
star is far from the black hole. In this regime, for a circular 
orbit with period less than P, the relationship between M 
and a can be derived from equation In cgs units, we 
have: 



M > 8 X 10" 



/ ar Y ( P_ 
UOi^cm/ 



(8) 



2.3.3 The apsidal motion bound 

The precession rate of peribothon for an equatorial orbit is: 



a^=(l 



(Ma^(l-e2))5 



(e.e:.. ICiufolini fc Wheeler! l|l995h 'l. 

Therefore, only for stellar orbits with large eccentricity 
can we observe the peribothon advance; at the same time the 
peribothon advance will give an estimate of the eccentricity. 
When the Cj is large enough to be observed, say uj ^ Comin, 
ignoring the small second term, we need: 



M > 



(1 - e^)LUi-niu 



(10) 



with a small correction if the black hole spins. Taking the 
peribothon precession period to be shorter than 10 years, 
and the eccentricity to be 0.5 for example, we have: 



M > 2 X 10 



2 / ar \- 
V10"cmy 



Mp 



(11) 



2.3.4 Lense-Thirring observability 



If the black hole is spinning and the orbital plane of the 
stellar trajectory is inclined with respect to the equatorial 
plane of the hole, the orbital plane will precess about the 
spin of the hole due to the dragging of inertial frames. This 
causes the prec ession of Q., the lo ngitude of the ascending 
node. We have l|Ciufolini fc Wheeler., 1995.) : 



9, ■ 



2SM^ 



a3(l-e2)! 



M > 



Therefore, in order to observe this Cl ^ fimin, 



2S 



(l-e^)3a: 



3/2 



(12) 



(13) 



The Lense-Thirring effect can reveal the spin of the black 
hole. 

For an orbit with e < 0.5, a Lense-Thirring precession 
period less than 10 yr, and S < 0.998, we will have: 



M>3xl0' 



/ ar \ ■ 
UO"cmy 



Ma 



(14) 
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Figure 2. Possible stellar trajectories around a slowly rotating 
black hole: the horizontal 'web' structure is the Keplerian accre- 
tion disc, ranging from 6 to 20 M. The four plots represent from 
top to bottom, stellar orbits with a. small eccentricity; b. large 
eccentricity; c. large distance from the black hole; d. close dis- 
tance to the black hole. Note that in each plot, the motion of the 
star is confined to one plane. The orbits are denoted by thick lines 
above the accretion disc plane, and thin lines below. 




Figure 3. Possible stellar trajectories around a rapidly rotating 
black hole with a = 0.998. The initial velocity and position of 
the star in each plot are the same as those of the corresponding 
stars' in Fig. |2] Note that the motions are no longer planar due 
to dragging of inertial frames by the Lense-Thirring effect. 
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2.3.5 Observability in the presence of gravitational 
radiation 

Gravitational radiation will carry off the orbital energy of 
the star. Eventually the star will be drawn in towards the 
black hole and plunge through the event horizon or be tidally 
disrupted. The gravitational radiation reaction time has to 
be reasonably long in order for the star to be observed. 
This time scale can be expressed as (iMisner et al.|[l973 ): 



for a non-spinning hole. 



(15) 



256 M^ALf(e) 

Here M^(^ M) denotes the mass of the orbiting star, and 
/(e) = (l+||e^ + §e*)(l-e2)-^ (|Peters fc Mathewdl 19631 ). 

For the stellar orbit to last for a useful finite time scale 
To, the M — ar relationship is: 



M ^at 



256roAf*/(e) 



(16) 



with a small correction if there is a nonzero S or in the 
strong field limit. 

Suppose we have observations of A'^agn AGNs and each 
one swallows TV* stars over its 10 Gyr lifetime. In order to 
find one example, we need: 

1 1 



TGR 



lOGt/r " iV* TVagn 



(17) 



Say N^, = 1O"^M*/M0 with M^, ^ IMq. If we observe 
10^", 10^, and 10* AGNs, the radiation reaction time is then 
respectively a thousand years, a million years, or a billion 
years, and the formula becomes: 

M ^ 4 X 10" ( I^Me, (18) 



M ^ 10' 
and 

M ^ 36 



4 / ar \ 
UOi^cmy 



( V 



M, 



©• 



(19) 



(20) 




10' 



10" 10" 10" 

stellar Orbit's Semi-major axis (cm) 



10^ 



Figure 4. Black solid lines (from right to left) - Line 1: To 
the left of this line, the period of a stellar orbit is less than 1 
year. Line 2: To the left of this line, the apsidal motion of the 
peribothon (taking arbitrary e = 0.5) has a period less than 10 
years. Line 3: To the left of this line, the Lense-Thirring effect 
becomes important, say, with a period less than 10 years. S = 
0.998 and e = 0.5. Green: To the left of this line, the star will 
plunge into the blackhole. 5 is taken as 0.998 to plot the ISCO. 
Dashed lines: Each line corresponds to a particular gravitational 
radiation reaction characteristic time scale. To the right of each 
line, this time scale is more than (from right to left)-Purple: 10® 
years; Yellow: 10^ years; Red: 10^ years. Dotted lines: To the 
right of these lines, the star will not be tidally disrupted. Each 
line (from right to left) corresponds to the tidal radius for: Blue: a 
sun-like star; Brown: a M8 dwarf star with radius 0.13-Rq, mass 
O.IM0; Dark Green: a white dwarf with radius O.OO9-R0, mass 
O.6M0. Example: The darkest shaded region corresponds to a 
stellar orbit for a sun-like star with orbital period less than 1 
year, apsidal motion's and Lense-Thirring effect's period less than 
10 years, and stable under gravitational radiation dissipation for 
more than 10^ years, taking e = 0.5 and S = 0.998. 



2.3.6 Tidal disruption bound 

The tidal radius is (e.g. lFrank fc Reed (|l976l )): 

- = 5xl0-f^f |^f^V*cm, (21) 

where i?* is the stellar radius. 

Therefore, for a star to avoid tidal disruption at a dis- 
tance Or to the black hole, the mass of the hole has to satisfy: 

«^(57w^i:)'t^"'«- '^^> 

For a main sequence star with mass less th an a solar 
mass would satisfy i?* ~ M°'* (|Kippenhahnlll993 ) . we have: 

which is: 

^^«>^^o'(Tofc^)'^« (24) 

for a sun-like star. 



^^4xl0^(-^) Mq (25) 
V lO^^cm/ 

for a dwarf star (M8) (taking R^, = O.13i?0, and M* = 
O.IMq), and 

^^^^>^i°''(io^)'^- (2^) 

for a white dwarf (taking i?* = Q.OO9i?0, and A/* = O.GMq). 
2. 3. 7 Summary 

We summarize these conditions in Fig. [J] in which the dark- 
est shaded region indicates a stellar orbit for a sun-like star 
with orbital period less than 1 year, apsidal motion and 
Lense-Thirring effects with periods less than 10 years, and 
gravitational radiation lifetime more than a million years. 

This plot demonstrates that, given our understand- 
ing of the distribution of stars in galactic nuclei (e.g. 
[Alexander fc Netzed[l99i : lGhez et al.l l^Ol')). it is quite rea- 
sonable to expect an example a sample of within 10* A'agn. 
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2.4 Inelastic collision of star and accretion disc 

In order to elucidate principles, we make a very simple model 
of a star-disc collision. Numerical simulations are needed to 
make a more sophisticated model. 

With the knowledge of the star and gas four-momenta 
at any space-time point, wo can calculate how much energy 
is radiated following the collision. We treat the star-disc col- 
lision {v-i,,Vgas ^ Vsound) as a totally inelastic collision of 
two massive particles, one with the mass of the star and the 
other calculated by multiplying the surface density of the 
disc by the collision area. Wc assume that the intensity of 
the flare is proportional to the energy deposited on the star. 
Here we assume that the stax goes through the accretion 
disc, 'cuts' a hole on that plane and carries away that paxt 
of the gas with it. The shape of the 'hole' is an ellipse of 
area A. 

Suppose that, before the collision, the four-momentum 
of the star is pj, and the four-momentum of the gas is pfg. 
After the collision, the four- momentum of the star is P2s, 
the four-momentum of the gas is pfgi and the total energy 
and average four- velocity of the photons are Eq and u". In 
the center of momentum frame of the photon, we have the 
conservation of energy-momentum of the system: 



pis + pig = P2s + P2g + Eau" , 

which gives, 

rusuis + niguig = msU2s + ragU2g 



(27) 



(28) 



Here rus is the mass of the star, and rUg is the mass of the gas 
which collides with the star. The w's are the corresponding 
four-velocities. 

Just after the collision, the star and gas move with the 
same four- velocity, so wf^ — U2g. In the frame of the star 
after collision, wfs and are (1, 0, 0, 0). We assume that 
the radiation is isotropic in this frame, so it™ also is (1, 0, 
0, 0) since isotropic photons will have no net mome To view 
a message, click on it.ntum. Squaring on both sides of the 
equation, we obtain 

2mgmsuiaU\ga. = —Eq — 2msEcC — 2mgEo — 2mgms. (29) 

We assume that the radiated energy and the rest mass 
of the gas are much smaller than the rest mass of the star, 
so we ignore terms like Eo^ and rrigEo compared to terms 
containing rUs- The remaining terms become: 



mgmsuisUiga — —rrisEo — mgnis. 
Therefore, the radiated energy is 
Eq ~ mg(-uisUiga - !)■ 



(30) 



(31) 



Here nig = a A is the mass of gas taken away by the star, 
where a is the surface density of the accretion disc, which is 
assumed not to change from collision to collision. Thus, the 
problem of calculating the radiated energy is transformed 
to calculating the size of the collision area in the disc rest 
frame. 

Here we assume that space-time is flat locally. The 
shape of the star in its rest frame is a sphere in 3-d space, 
ignoring the effect of gravity exerted by the central black 
hole and rotation. In flat space-time, the major and minor 
axes of this ellipse are simply i?* and R*/cos{6g), where 9g 



is the angle between the normal vector of the plane and the 
velocity of the star in the rest frame of the gas in the disc. 

Define the frame of the gas such that the accretion disc 
is in the xy-plane, and the normal vector is the z-axis. Sup- 
pose that a and b arc the major and minor axes of the impact 
ellipse. In the gas's initial rest frame, its area is: 

A=^ (32) 

cos Og 

In the Boyer-Lindquist frame, there is a Lorentz con- 
traction along the direction of the motion of the gas in the 
plane of the disc, so the area becomes 

A=^^ (33) 

cos 9 gjg 

where 7g is the Lorentz factor of the gas. 

Clarifying the notation: in the gas frame, the four- 
velocities of the gas and the star are 



(1,0,0,0), 

/ 1 X y z \ 
'Jsgyl, Vsg, '^sgi '^sgj- 



(34) 
(35) 

Here the Xab means object a's quantity x, in b's frame, and 
from these two equations we have 



^g Usa , 



COS^g = \Vsg\/Vsg, 



7s9 



(36) 
(37) 

(38) 



In the Boyer-Lindquist frame, the disc lies in the r — </> 
plane, and the four velocities and the unit normal vector of 
the disc plane are 



Ug = ^g(l,Vg,Vg,v"g), 

a. / 1 X z y \ 

n" = S"^ (0,0, 1,0). 

Using special relativity, we have 

z 



79(1 



(39) 
(40) 
(41) 

(42) 



product. Vg is the 6 component of Vs in the Boyer-Lindquist 
frame, so 



Also, the inner product of Ug and Us gives 

Ug Usa 



1- Vg ■ Vs 

Therefore, 
cos Og = - 



(M^Msa)-\/l - 7s/ 



(43) 



(44) 



(45) 



The collision area calculated in terms of four-vector 
scalar products is then: 



^ {UgUsa) 



yielding the total collision energy as 

2 Ug Uscx ^" 1 



-a'nRi,- 



^g\Usna\ 



-\J {UgUs 



- 1 . 



(46) 



(47) 
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2.5 Ray-tracing 

The energies of tiie pliotons in tlie o bserver's frame are 
sliifted by a redsliift factor g — E/Eq (| Cunningjiaml [l975l ) 
compared to tliose in the emission frame. E is the energy 
in the observer's frame and Eq in the star's inertial rest 
frame. Noting that E = —p"ua and Eq = —pt, we have 
g = pt/{p"ua), where Ua is the four- velocity of the star. Us- 
ing conservation of phase-space, it is possible to show that 
this causes an ex tra inte nsity shift proportional to the fourth 
power of energv. fLindau ist (1966) shows that bolometric in- 
tensity / = g'^Io- For spectral intensity, I^, — g^I^o at the 
transformed frequency. 

We ray trace the photons backwards in time from the 
observer to the disc. Each photon has a different propaga- 
tion time to the observer, due to their different emission sites 
on the disc, their initial momentum, and the Shapiro time 
delay. The different time delays of the photons contribute to 
the quasi-periodity of the flares. By collating the mapping 
between the sky-plane and the surface of the disc, it is pos- 
sible to obtain the light delay time and redshift factor for 
each image pixel. Using this information, one can calculate 
the flare pattern for any stellar trajectory. This is efficiently 
done by sorting the pixels by light delay time and scanning 
only those pixels known to be within an impact region. 

We use a top-hat function for the shape of the flare in 
time. By using such a simple functional form, we seek to 
describe the timing accurately. However, this would have to 
be convolved with a more realistic flare profile to produce a 
light curve to fit observational results. 



3 RESULTS: ARRIVAL TIME OF THE FLARES 

With the knowledge of the stellar trajectory, the collision 
model, and the ray trajectories, we calculated the arrival 
time profile using the mapping between flare locations and 
the observer. The signals exhibit quasi-periodicity for two 
reasons. Firstly, due to the precession of the stellar orbit, 
the time separations between two consecutive collisions are 
different. Secondly, because the collisions happen at different 
locations each time, the light travel times for the signals 
differ. The nature of this quasi-periodicity depends on the 
mass and spin of the SMBH and the orbital elements of the 
stellar orbit, as summarized in Table [1] Some example plots 
are shown in Figs. [5]and|6l 

Eight parameters excluing the epoch To are needed to 
characterize the time series in general: the mass of the black 
hole, the spin of the black hole, the observer's inclination 
angle w.r.t. the spin, and the orbital elements of the the star. 
In addition, two of the angle parameters are degenerate due 
to the symmetry of the model to rotations about the axis in 
the direction to the observer. 

Finally, in highly symmetric orbits, other parameters 
may be degenerate. The eccentricity is related to the apsi- 
dal motion, and can be determined from the changing rate 
of the period. These two obviously become degenerate for 
circular orbits. The spin and the inclination are coupled by 
the Lense-Thirring effect, and if the black hole does not spin, 
the orbital inclination is not measurable. 

If the peribothon is midway between two hits, the en- 
ergy released in successive collisions will be similar. However, 



Table 1. Behaviors of quasi-periodicity according to tlie param- 
eters 



Parameter Effects on quasi-periodicity 

Or The larger the semi-major axis of the stellar 

orbit, the more periodic the signals are. 

e A large e causes the interval between 

two consecutive events to differ. 

i Orbits with higher inclination are 

affected more by the black hole spin. 

f2 Different Roemer time delay 

and Shapiro time delay 

LO Ratio of the height of flares 

6 Large observation inclination angle 

increases flux intensity variability. 

M Changes overall time scale 

S Causes the time separation 

between signals to vary 



due to our assumption of isotropic emission in the star's 
frame, if the star is approaching the observer in a collision, 
the observed collision will be brighter due to Doppler boost- 
ing. These effects mean that ffares will vary in brightness. 
Unfortunately, our model of flare flux is quite primitive, but 
this does not affect the use of the flare timing in determining 
system parameters as long as the flares can be individually 
distinguished. 

Finally, note that the small spikes at the base of the 
large spikes most noticeable in Fig. [5^ are flux from higher- 
order images due to photons partially orbiting the black hole 
whilst propagating to the observer. This can provide extra 
information about the inner extent of the accretion disc due 
to it possibly blocking this flux. 

3.1 Dominant parameters 

ar,M,e, and Q. basically determine the temporal behavior 
of the signals while the other parameters give small correc- 
tions to the quasi-periodicity or change the flux variabil- 
ity. Therefore, we investigate these three parameters in a 
Schwarzschild case since the spin is not that important here. 

3.1.1 Size and eccentricity 

Figs. [5]d and [5}; show the observed flares resulting from 
stellar orbits with smaller and larger respectively, corre- 
sponding to orbits in Figs.[2}5 and [2};. As can be seen, when 
the size of the star's orbit becomes large compared to the 
gravitational radius Rg, the signals becomes more periodic. 
This is because when the size of the star's orbit is large, or 
on average the star is far away from the black hole, general 
relativity plays a rather unimportant role. The star basically 
moves as on a Newtonian orbit, giving periodic patterns of 
the signals as the precessions is small. Note that the mass 
of the black hole only appears as a scale in the problem, via 
the size of Rg. 
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Figure 5. These are the fluxes from the flares corresponding to 
the steUar trajectories in Fig. [2] The x-axis is time in the units 
of The y-axis is the fluxes of the flares in arbitrary units. 

From top to bottom, in [2^, the orbit is basically circular, and 
consecutive flares are separated by similar time intervals. In [2b. 
the orbit is more eccentric, and flares are unequally spaced. The 
flares from[2j; are much more periodic than from[2ji due to the 
star spending more time away from the black hole. The bottom 
case shows how the secular effects of precession can add up to 
large changes in the flare pattern. 



Figure 6. These are the fluxes of the flares in Fig. [3] The spin of 
the black hole affects the orbit via inertial frame dragging. This 
introduces extra precession which can be seen as modulation of 
the flaring pattern. In principle, this allows the determination of 
the black hole spin from observations of star-disc collisions. 
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Figs. [5^ and[5]3 show the flares from the two trajectories 
in Figs. [2^ (low e) and[2l3 (high e). 

Consecutive signals may be grouped in multiple peaks. 
The class of trajectory that does this is described in Fig. [T] 
As can be seen, the time intervals between a group of three 
consecutive peaks can be much shorter than that between 
the third and fourth signals, and this pattern will repeat 
itself. This gives a four-peak feature, as shown in Fig. [SJj. 
Note that if the star goes between the hole's event horizon 
and the inner radius of the accretion disc, or beyond the 
outer radius of the accretion disc (as in Fig. [5};), we will 
miss some hits and may observe only double peaks. 



3.1.2 Longitude of the ascending node 

To examine the role that fl plays, we use the same stellar 
trajectory (as in Fig. [^Jj) and rotate it about the z-axis. For 
the original orientation, the hits' longitudes are about and 
TT, namely, in front of and behind the black hole as seen by 
an observer. Then for a second orientation which is a tt/2 
rotation of the orginal one, the star goes through the disc 
on the sides of the black hole with = n/2 and 3tt/2. And 
for a third orientation which is a 7r/4 rotation, the situation 
is half-way between these extremes. 

The corresponding fluxes with a 5° (top) and 45° (side) 
observing angle are plotted in Figs. |8] and [O] In the first fig- 
ure, we find the time interval of any double peak is nearly 
independent of its set-up. However, in the second figure we 
see that the intervals between peaks are altered. This is be- 
cause when we are observing from above, the longitude of 
the ascending node does not affect the propagation time of 
the signals due to circular symmetry. However, when we are 
observing from the side, the longitude is an important factor 
in determining the time delays. The flares behind the black 
hole need more propagation time than those in front - the 
Roemer delay. The order of this time correction is compa- 
rable with the time correction due to the precession of the 
orbit, thus causing the noticeable change in the duration 
of each double peak. Also, the Shapiro delay times are dif- 
ferent, which is caused by the light propagation in curved 
space-time taking a slightly different path than in flat space- 
time. Both time delays are the same, however, for signals 
from the sides due to symmetry since it is non-spinning. 



3.2 Orbital inclination i and spin of the black hole 

S 

The inclination of the stellar orbit is relevant only when the 
spin of the black hole is large. For a Schwarzschild black 
hole, the inclination plays no role in determining the quasi- 
periodicity or the flux of the signals. Also, the effect of spin 
of the black hole is not prominent unless the orbit is highly 
inclined. 

The spin makes the orbit less periodic via the Lense- 
Thirring effect, and it no longer is conflned to a plane. How- 
ever, the change to the flare time-series is relatively simple. 
The result is similar to that from a non-spinning black hole, 
with the frame dragging inducing an extra modulation of 
the flare fluxes and timing. This can be seen by comparing 
the fluxes in Figs. [5]and[6l 




Figure 7. This figure shows how multiple flares can be produced. 
The star's trajectory can be strongly affected by the black hole's 
gravity. This can cause more than two crossings of the accretion 
disc plane per close approach. Depending on the inclination of 
the orbit relative to the accretion disc, the extent of the disc, 
and the amount of doppler boosting, a varying number of closely 
separated flares may be seen; i.e., for the case where the flares 
are produced at locations 1, 2, 3, 4;, there may be two bright, and 
two dim flares per close approach. For the cases 5, 6, 7; there 
may be either two bright and one dim, or two dim and one bright 
flares. Also note that extreme gravitational lensing may produce 
higher-order images of flares. This may introduce more features 
in the light curve if the orientiation is such that the accretion disc 
does not obscure them. 

3.3 Observer's inclination angle 6 and argument 
of the peribothon ui 

These parameters only marginally affect the temporal be- 
haviour via differences in the Shapiro and Einstein time de- 
lay. However, they can greatly affect the flux intensity. Figs. 
1101 and [TT] show how the strengths of fluxes change as the 
observer's inclination angle changes. When the inclination is 
large, some signals are greatly magnified via Doppler boost- 
ing, depending on their locations on the accretion disc. 

The total energy emitted by the flare depends on the 
velocity of the star in the gas frame and also upon the rela- 
tive direction to the line of sight. If one assumes models of 
accretion disc structure and flare physics, this provides extra 
information either to determine the parameters of the disc, 
or to constrain the orbit of the star since Doppler boosting 
is a sensitive function of angle. 

Also note that if two have similar parameters except for 
the argument of peribothon, flares closer to the peribothon 
will be more intense intrinsically. This is because the star 
is moving at a higher speed there. This results in larger 
variability in the flux intensity proflle. 
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Figure 8. These fluxes are for the orbits in Fig. ??, when the 
observer is observing almost from top down at a 5° angle with the 
spin axis of the black hole. The time intervals of any double-peak 
is roughly the same in any figure. The spin of the blax;k hole is 
taken to be for simplicity in this figure and the next one. 



Figure 9. The associated orbits are the same as in the previous 
figure, but the observing angle is changed to 45°. Light travel 
time delay makes the observed time between flares uneven when 
the flares come from the front and back of the black hole. Also, 
the effects of Doppler boosting are enhanced. 
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Figure 10. This figure shows how the relative strength of the 
signals change with the observer's inclination angle. We used the 
star's orbit in Fig. [2J3 with a Schwarzschild black hole in the 
center. The five figures correspond to observer's inclination angles 
of 5° (from the top), 30°, 45°, 60°, and 85° (from the side). As 
shown in the figures, when observing from the top, the scales of 
the fiuxes are relatively similar, while as the inclination angles of 
the observer increases, some fluxes' scales start to be increased 
or decreased. When we are observing from the side, some signals 
are greatly magnified. 



Figure 11. This figure shows the same as Fig. 1101 but using 
[3I3 and a rapidly spinning black hole. The higher the inclination, 
the larger the maximal Doppler boosting effect for the correctly 
aligned flares. 
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4 FROM THE SIGNALS TO THE 
PARAMETERS 

The parameters we would like to measure from observa- 
tions of star-disc collisions are M, S,ar,e,Q,LL!,i,To, and 
6. However, as we discussed in the previous section, un- 
der different circumstances, some of these observables are 
degenerate: they are coupled and cannot be determined in- 
dividually. Here we study analytically which parameters or 
combinations of the parameters could be determined from 
the observations of quasi-periodicity of the signals for the 
Schwarzschild black hole case. 



4.1 Slowly spinning hole 

In this case we can d escribe the space-time 
the Schwarzschild metric. iBlandford fc Teukolsb " 



Damour fc Deruelld lll985h . iDamour fc DeruelL 



using 



197e), 



19861 ) ■ 



Damour fc Tavlod (ligod ). and lDamourl (119931 ) explicitlv caL 



culated the post-Newtonian timing formula of binary sys- 
tems. We employ their formulae since our binary system 
composed of a Schwarzschild black hole and an orbiting star 
also behave as a binary system in the moderately strong- 
field regime of relativistic gravity. They calculated the time 
delays of the arrival time of signal as; 



D-Ta=T + AR{T) + AEiT) + AsiT) + ^-j Pj (48) 

up to an arbitrary additive constant, with T as the intrinsic 
emission time, and each physical quantity is evaluated at T. 
(1) Roemer Time Delay; 



Ar{T) = -fc-r*. 



(49) 



where fv is the coordinate postion vector from the star to 
the barycenter of the binary system, which is the SMBH 
since we assume its mass is much bigger than that of the 
star, and k is the unit vector opposite to the light of sight. 
(2) Shapiro Time Delay; 



As(T) = -2Mlog [k ■ (rV - rBu) + \r, - rBni 



(50) 



where vbh is the position vector from the black hole to the 
barycenter of the binary system, which, in this case, is 0. 
(3) Einstein Time Delay: 

Ai3(T) = 7sinu. (51) 

Here we have 7 — 2Me/nar, with the mean motion n — 
^ M/af{l — 9Af /2ar), and the eccentric anomaly u satisfies 
n{t — To) = u ~ e sin II, with t as the emission time of the 
signal. 

The star hits accretion disc when; 



k ■ fi, 



r* cos a — ±r* cos 6 cos Q, 



(52) 



where a is the angle between the line of sight and f*. The 
sign is different for the ascending or descending node. 

At the ascending node, rv makes an angle u with the 
peribothon. Therefore we have: 



(53) 



where eg is a constant with the value e(l 4- 2M/ar), and 
A^g (u) is defined as; 



Aag{u) :— 2arctan 



'n + eg-^ 




tan — 


\1 — egy 


' 2 



with the property; 



cos {Aeg{u)) = J 



COS u ~ eg 



I COS u 

Therefore after calculation; 



COSM 



and 



ar(l - 



el) 



l + egcos{jf^)' 

Similarly, at the descending node, 

^ ^ "r(l~ei) 

l_e,cos(^)' 

Thus, we obtain: 

ar{l - el) 



k ■ f 



sin {Aeg{u)) 
we obtain: 



1 ± Br COS ( 

th 

so, 

iinuy/l — eg 
1 — eg cos u 



cos 6 cos f2, 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



where — stands for the ascending, and + stands for the 
descending node. Also, using 



l±eg cos(Yf-) 
Therefore, the time delays are 
ar{l 



Ar{T) = T 

As(T) : 

and 
Ab(T) 



el) 



l±ee cos ( j-^) 
ar(l 



cos 9 cos 57, 



-2Mlog 



2Me 



- eg) 



1 ± ee cos I 



(1 ± cos cos Q) 



na,- l±eg cos ( 



1+LU 



(60) 
(61) 

(62) 
,(63) 

(64) 



Define a as the angle between the line connecting the 
SMBH to the observer and the line connecting the SMBH 
to the ascending node, then cos 6 cos Q. = cos a. We cannot 
decouple 6 and Q., since the system is rotationally invariant 
about the line at which the stellar orbital plane intersects 
the accretion disc. In addition, define u = tj/(l -I- uj). u 
equals the initial argument of peribothon uj if there is no 
precession. Then the three time delays become; 



Afl(T) = 
As{T) = 

Ae{T) = 



l-e\l + q)' 
-r-. — TT^ — ^ — cos a 

1 ± e(l -f q) coscij) 

r.(l-e^(l + g)2) 



-a^ijlog 



1 ± e(l + g) cos UJ 



(1 ± cos a) 



2^ ^l-e2(l + g)2 



[l-2^^l± e(l -I- g) cosii 



(65) 
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where q = 2M/ar is dimensionless quantity, which compares 
the scale of the black hole and the scale of the stellar orbit. 
In addition, the precession rate of peribothon and the orbital 
frequency become 



4ar(l - e2) 

2 



(66) 



These define the observables obtainable from pure timing 
measurements. 

Theoretically, one would only require six flares to solve 
for the reasonable parameters of the system: T*, g, e, cos a, 
and uj. There are five unknowns, and five equations via (|65|) 
and (|66|l for them, together with an extra equation/unknown 
for the time initialization Tq. These theoretically can be 
measured and simultaneously solved to determine the sys- 
tem. 

Unfortunately, the above equations cannot be directly 
used to calculate the parameters of the system from obser- 
vations. They are only approximations, and do not include 
the effects of the black hole spin on the light delay time. 
In addition, they assume that the orbit is an ellipse which 
is incorrect close to the black hole. Note that the orbital 
trajectory can be integrated a nalytically in terms of elliptic 
functions l|Fanton et al.|[l997l ). However, these expressions 
cannot be easily inverted. Thus, a solution will require some 
form of numerical fitting of the observations and more than 
the minimal number of flares as data points. 

The best approval appears to be hierarchical fitting. 
Some parameters are much easier to determine than others. 
By looking at the modulation of the fiare signals, it should 
be possible to obtain the period and precession rates. With 
these, it seems to be possible to reduce the size of the multi- 
dimensional space enough that a direct search of possibilities 
for the remaining parameters is possible. 



4.2 Rapidly spining holes 

In this case, we must use the full Kerr metric, and only 
limited progress can be made analytically. An analysis of 
observational data would undoubtably involve numerical in- 
tegration of the orbit and ray equations. 

When we add the spin to the group of the parameters, 
one more equation needs to be added to solve the problem in 
principle, which is the Lense-Thirring precession rate. The 
precession rate of orbital frequency is: 



2r*(l--e2)f 



(67) 



The other five equations need to have corrections due to the 
spin too, but we can express two of them analytically: 

1 / ■iqy/2q q^S 



±1 



(68) 



while the other three time-delay formulae are extremely 
hard to calculate analytically. Fortunately, we have another 



way to estimate the spin parameter besides using numerical 
methods. 

As we have shown in Fig. [6l the spin of the black hole 
would transfer power to higher harmonics in the time se- 
ries, which increases the difficulty of studying the spin of 
the hole through the timing analysis. However, due to the 
precession of the longitude of the ascending node and de- 
scending node, fiares triggered at different locations of the 
disc would have different Doppler shifts in their intensity. 
Since the peribothon advance affects the intensity change 
marginally, we have a reason to believe that the flare inten- 
sity modulation is basically caused by the Lense-Thirring 
effect and the period of this modulation is related to the 
spin of th e black hole. 

From lCiufohni fc Wheeleil (|l995l ). we have: 



S = 



2M2 



17. 



(69) 



If we define the Lense-Thirring precession period to be Pn = 

2-k/Q,, 



S = 



■Kar{l - e^)i 
M^Pn 



(70) 



Tests using the modulation period of the flux in Fig. 
[H]and the size and eccentricity of the stellar orbit from the 
simulated orbits in Fig. [3] have been performed in which the 
spins are all calculated to be correct in order of magnitude. 
Therefore, the period of the modulation of the flare fluxes, 
together with the size and shape of the stellar orbit, and 
the mass of the black hole inferred by the orbital period, 
can provide a probe into the spin of the hole, though more 
studies need to be done on how the spin can be calculated 
more accurately directly. 



5 OBSERVATIONS 

5.1 X-ray flares from tidal disruption of stars 
near slowly accreting SMBHs 

Stars on plunging orbits in a galactic core m ay approach 
close enough to the SMBH to lose mass (jRee^ Il988l : 
lEvans fc Kochaneklil989 ). A fraction of the outflowing gas 
may be captured and accreted by the black hole, producing 
out bursts of radiation . 

iKomossa fc Bade! (|l999l ) reported non-recurrent X-ray 
flares from NGC 5905 and IC 3599. The former is an inactive 
galaxy, and a tidal disruption event with peak luminosity 
exceeding Lx ~ 10*^ e rg s~^ flts the da t a best . Using the 
Eddington luminosity, iKomossa fc Bad"3 (|l999l ) calculated 
that the lower bound of the SMBH mass would be 10^ M©. 
This outburst lasted for several months, consistent with the 
theoretical cal culation of iReei l|l990l ). 

Moreover, iKomossal 1 200 ll ) reported a few more X-ray 
flares from inactive galaxies, all of which ould be interpreted 
as tidal disruption events. These soft X-ray flares have large 
luminosities, high degrees of var iability, and the abs e nce o f 
optical signs of Seyfert activity. IKomossa fc Merritd (|2008D 
further discussed this scenario in the environment of recoil- 
ing SMBHs. 
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5.2 X-ray Periodicity from the AGN in galaxy RE 
J1034-I-396 

iGierlinski et al.l (120081 ) reported the observed periodic 
XMM-Newton X-ray signals from a nearby active Seyfert I 
galax y RE J1034+396. As shown in Fig. 1 in.Gierhnski et afl 
l|20od ) ■ the signals have a periodicity around 1 hour (more 
precisely, at 3730 ± 130 s). The X-ray flux has a modula- 
tion of around 10% and is identified with the nucleus on 
bolometric grounds. 

The data were divided into two segments. For the first 
seven periods, there are phase shifts and extra maxima. 
Then after to = 25 ks, in the second segment , there are 
about 16 coherent cycles. IGierlinski et al.1 (|2008D folded the 
second segment and showed that the X-ray modulation could 
be well fitted by a cosine function. 

We consider two scenarios to explain the observed X-ray 
periodicity and sinusoidal light curve. The first is that the 
star orbits equatorially inside a thick disc, and emits X-rays 
due to its velocity relative to the disc. The relative velocity 
might result from the fact that the disc will be pressure- 
supported and sub-Keplerian, the gas within the disc can 
have an inward velocity, or the star orbits with a nontrivial 
eccentricity. The rapid inflow would preclude the creation of 
a gap in the disc. Another scenario is that instead of having 
photons produced directly by the star-disc interaction, we 
can have a star slowly overflowing its Roche lobe as it orbits. 
This would be a source of accreting gas and the modulated 
emission may be produced at a variable hotspot where the 
accretion stream from the Lagrange point hits the disc. If 
this is the case, predictions of the the evolution of the stellar 
orbit can be made. Further discussion of this model and RE 
J1034-f 396 will be described elsewhere. 



5.3 Periodic optical outbursts of OJ287 

The BL Lac object OJ 287' s optical outbursts exhibit peri- 
odicity of roughly 12 years (iNilsson e t al. (2 00^, and refer - 
ences therein), as shown in Fig. 1 in iNilsson et al.l (|2006l ). 
The last two outbursts also show a double-peak feature of 
roughly a year separation. 

Many models have been suggested, and one of them pro- 
posed by Lehto & Valtoncn (1996) suggests that these peri- 
odic flares can be fitted by the orbit of a secondary 10* 
SMBH around a primary 1.7 x 10^° Mq SMBH. The orbital 
period is 12.07 years, and two flares are produced when the 
secondary SMBH crosses the accretion disc of the primary 
one due to thermal bremsstrahlung. As shown in Fig. 2 in 
INilsson et al.l l|2006l) . the time separation between flare 1 
and 2 is roughly a year. This model is equivalent to ours 
with a 1.7 X 10^°Mq Schwarzschild central black hole. Since 
the secondary black hole is relatively far from the center, 
the orbit is nearly periodic. Here the orbit has a precession 
of about 33° per orbit. More observations will be necessary 
to see if this model is highly predictive. 



5.4 Infrared flares from the Galactic center 

Our Galactic center hosts a supermassive black hole of mas s 
4.1 ~ 4.5x 10^ MfTi llGenzel et al.l (|2003l ): lGhez et all ||2008D : 
iGillessen et al.l (|2009l ) : and references therein) , which gives a 



gravitational radius ijg = M ~ 18 light seconds. The inner- 
most stable orbit, which ranges from 1.24 to 6 Rg depending 
on the unkown spin of the black hole, will be about 0.1 ~ 1 
light min ute. The radio size of the black hole , as detected 
by VLBI (|Bower et al.ll2004l : IShen et al.ll2005l ). is less than 
10 light minutes. 

Two groups mappe d the central sta r s in t he Galactic 
center, as in Fig. 2 in iGenzel fc KarasI |20o3). The clos- 
est observ ed stellar orbit peribothon is about 17 light hours 
in radius (ISchodel et al.ll2002l . l2003l : iGhez et al.|[2003l . l2005l : 
lEisenhauer et alT2005^ . These observable stars are bright 
OB stars. For a given supermassive black hole, the tidal 
radius Rt of an orbiting star is proportional to i?*/Aft^'^^ 
(|Reeslll988l ). In the l ower mass range of main sequence stars, 
we have 7?* oc M^'^ l|Kippenhahnll 19941 ). Therefore, for main 
sequence stars Rt oc Mf ''^. This means the less massive a 
star is, the denser it is, and the closer it can get to the cen- 
tral black hole. Therefore, a sun-like star can get as close as 
around 5 light minutes to the center of the black hole in the 
Galactic center; a dwarf star can get to around lORg; white 
dwarfs and neutron stars have higher densities than main 
seq uence stars, and will cross the horizon intact. 

iGenzel et al] (|2003l ) detected quasi-periodic infrared 
flares using ESO-VLT, as shown in their Fig.2. The signals 
have a 17 min period between consecutive peaks and a 1 day 
period between groups of peaks. If we try to explain these 
flares with our star-disc collision model, we need to invoke 
a highly eccentric stellar orbit around a Kerr SMBH with 
an orbital period ~ 1 day. This period would allow almost 
all stars to escape from tidal disruption. However, we do not 
have any observation of stars that close to the hole, and there 
are 3-5 consecutive peaks in each group, which we cannot 
explain with our model. Suppose a hot spot is produced in 
each collision and it orbits around the hole for several cycles. 
The consecutive peaks should then follow a power-law curve 
instead of the rather symmetric distribution in Fig. 2e in 
Genzel's paper. Therefore, a star-disc collision model seems 
an unpromising explanation of this quasi-periodicity. 



6 DISCUSSION 

In this paper, we have explored observable effects that can 
arise when a normal star follows a bound orbit close to the 
event horizon of a massive black hole in a galactic nucleus. 
We have paid special attention to the quasi-periodic signal 
that may be detected in X-ray observations if the star cre- 
ates an X-ray pulse every time it crosses the disc plane. We 
have also considered, briefly, alternative forms of signal mod- 
ulation, especially associated with more subtle star-disc in- 
terac tions involving hot spots and tidal shocks (|Spruit et al.l 
[l983)- 

If the orbit is relativistic, the Doppler effects will al- 
most guarantee a measurably modulated signal. 

Although the detection of effects like this would be of 
considerable interest in its own right, it would open up two 
more far-reaching lines of inquiry. The first of those concerns 
the testing of general relativity in the strong field regime. To 
date, all quantitative tests of relativity have been executed in 
the weak field regime. (We exclude interpretations of obser- 
vati ons such as a symmetrically broad iron fiuorescence lines, 
e.g.lMilleil ll20q7^ and QPOs in stellar b ack hole discs, e.g. 
IWagonej (|l999l ): IWagoner et al.l l|200ll ).') We have described 
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the types of time series that could be measured and shown 
how these can be interpreted initially as measurements of 
orbital and black hole parameters and then, given sufHcient 
high-quality data, as prescriptive tests of the general rela- 
tivistic calculation of orbits in a Kerr space-time. Such a 
test would represent a major advance, though it would not 
be as powerful as the measurement of a gravitational radiar 
tive waveform from a major black hole merger. 

The second line of inquiry is that observations would 
allow us to quantify the rate of extreme mass-ratio inspi- 
rals (EMRI) in galactic nuclei. (The discovery of 100 short- 
lived OB stars in the Galactic center underscores how little 
wo understood about the co-cxistcnce of stars and massive 
black holes in galactic rmclci.) This would, in turn, relate 
to understanding the role of stars in the behavior of mas- 
sive accretion discs. It would also relate most directly to the 
planning for observations such as LISA that are designed to 
detect gravitational signals from stellar inspirals. 

Of course, most exciting of all would be to understand a 
stellar orbit well enough that the stars "quietus" - either at 
the ISCO or through rapid tidal destruction - could be pre- 
dicted and observed. Given our current ignorance, we cannot 
rule out such an opportunity arising within our lifetime. 
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